Kagome Lattice Antiferromagnet Stripped to Its Basics 
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We study a model of a spin S = 1/2 Heisenberg antiferromagnet on a one dimensional lattice with 
the local symmetry of the two dimensional kagome lattice. Using three complementary approaches, 
it is shown that the low energy spectrum can be described by two critical Ising models with different 
velocities. One of these velocities is small, leading to a strongly localized Majorana fermion. These 
excitations are singlet ones whereas the triplet sector has a spectral gap. 
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The famous kagome lattice antiferromagnet still re- 
mains largely a mystery after a decade of extensive stud- 
ies. This system exhibits both frustration and low coordi- 
nance and classically it has infinite continuous degenera- 
cies. Local distortions allow it to explore its many ground 
states with no cost in energy and lead to a very specific 
linear spin wave spectrum with a whole branch of zero 
energy excitations |Q. In the quantum case (S=l/2), 
the system is likely to be a spin liquid with a gap for 
the magnetic excitations 0. For finite samples, the sys- 
tem has a huge number of singlet states below the first 
triplet H which is a rather unexpected feature for a two- 
dimensional quantum antiferromagnet. Moreover, the 
analysis of the specific heat shows the existence of un- 
usual low- lying excitations Q. The presence of these 
low-lying singlet excitations below the spin gap gives a 
picture of an intriguing spin liquid that deserves under- 
standing from a general point of view. 

In this paper, we shall study a one-dimensional model 
of spin S = 1/2 Heisenberg antiferromagnet on a lattice 
presented on Fig. ^ which retains the local symmetry of 
the kagome net. Insight into the behavior of the two- 
dimensional kagome antiferromagnet might be gained by 
investigating this simplified system. In particular, the 
one dimensional model may help us to identify the slow 
degrees of freedom of the problem. 

According to the conventional wisdom the model we 
study may be viewed as a version of a three-leg spin lad- 
der which in the low-energy limit one should expect to 
fall into the universality class of the S = 1/2 Heisen- 
berg spin chain. The latter means that the low-lying 
excitations (often called spinons) are represented by one 
gapless bosonic mode (in the language of the theory of 
critical phenomena this means that the central charge of 
the model is equal to C = 1 ||). However, the frustra- 
tion may play its tricks and, as we shall see, a somewhat 
different scenario is realized. Namely, the bosonic mode 
decouples into two modes of real (Majorana) fermions 
(C= 1/2 each) having different spectra. 
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FIG. 1. One dimensional version of the kagome lattice. 

The Hamiltonian for the model shown in Fig. |l| may 
be written 

7~L = J]] / (Sa,2j ' S Qi 2j + l + S Q: 2j + l • S a: 2j+2) 

0=1,2 j 

+ J±S 2j+ i ■ (S a ,2j + S a ,2j+l)] (1) 

where s 2j+ i are the "middle spins", and S aj are the 
"rail spins" , the index a taking two values [a = 1,2), one 
for each rail. All interactions are antiferromagnetic and 
we shall also consider the case where interaction between 
spins belonging to the rails of the lattice is much stronger 
than the interactions with the middle spins: Jii >> J±. 
Under this condition, one can describe the chains in the 
continuous approximation, representing the spins on each 
chain as a sum of ferromagnetic (M a ) and staggered (n a ) 
parts: S a ,j — > M a (x) + (— l) J n a (a;). It is crucial for our 
analysis that the middle spins interact only with the fer- 
romagnetic part of the magnetization. The latter has the 
following remarkable property - it can be written as a sum 
of "currents" M a = J a + J a , where the currents satisfy 
the same commutation relations as bilinear combinations 
of left- and right-moving Dirac fermions m : 



1 



— Ra, a <r a pRa,p/'2,'fa. = L+ a <J a pL a ^ /2 



(2) 



where a are the Pauli matrices. The corresponding alge- 
bra (SU(2)i) is called the level k = 1 SU(2) Kac-Moody 
(KM) algebra where k refers to the number of species of 
spin- 1/2 fermions (one in the given case). Now we no- 
tice that the interaction includes the sum of two SU(2)i 
spin currents J = Ji + J 2 which, by definition, is a level 
k = 2 current. These currents have the same commu- 
tation relations as currents of three Majorana fermions 
(the explicit expression is given later) QJ^]. Since each 
Dirac fermion can be represented as a linear combination 
of two Majorana (real) fermions, it is clear that only 3/2 
of the low-energy degrees of freedom of rails are involved 
in the interaction. Obviously the other decoupled degrees 
of freedom remain critical. Using the results of Ref. ||] 
we can represent the continuous limit of the Hamiltonian 
(El) as follows: H = H S +H t with 
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H s = - — dx (r d x r - l dj ) 
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dx [3 2 (x) + 3 2 {x)} 



Y^9 J dxS(x~ aj) Sj ■ [J(x) + 3{x)}. 



(3) 
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Here v ~ J\\ is the spin velocity which we shall later put 
to unity and the coupling constant is g ~ Jj_. The Ma- 
jorana fermion (ro,io) describes non-magnetic, gapless 
excitation. It represents Ising degrees of freedom that 
do not interact with the central spins and are associated 
with a discrete Z 2 interchange symmetry between the 
surface chains. We shall study the Hamiltonian (Q) using 
several alternative approaches. First, we shall study an 
integrable deformation of model (Q) . To make sure that 
the exact solution reproduces the qualitative features of 
the spectrum, we shall consider an anisotropic version 
of (||) bearing in mind that it may simplify in the limit 
of strong anisotropy (the so-called Emery-Kivelson limit 
|1). Finally, a direct mean field treatment of the lattice 
in the isotropic limit will be performed. 

Exact solution of a related model. Since the cur- 
rents satisfying the level k = 2 SU(2) KM algebra can 
be represented as fermionic bilinears, the effective the- 
ory for the central spins corresponds to a special version 
of the two-channel (k — 2!) Kondo-lattice model where 
"electrons" do not experience backscattering. A similar 
model (with k = 1) was considered in Ref. [M; the result 
is that integrability is achieved if one adds an additional 
interaction: 
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Ht + g / dx 3(x) ■ 3(x) 



(5) 



where the ratio g' / g is fixed. The Bethe ansatz equations 
for model (J5J) are derived in the same way as in Ref. |lC| ] 
with the only difference that "conduction electrons" now 
belong to the spin-1 representation of the SU(2) group. 
The result is 



M 

[e 2 (x a + \/g)e 2 {x a - 1/ g)] N ef° (x a ) = e 2 (x a - x b ) 

6=1 

(6) 

where 27V and Nq are the numbers of spins on the rails 
and in the middle, respectively, M is the number of up 
spins (M = 2N + N /2-S z ) and e n (x) = (x-in/2)/(x + 
in/2). The total energy takes the form 

E = I He2(x a + l/g)/e 2 (x a - 1/g)}. (7) 

a 

Following the standard procedures of the Bethe ansatz, 
we derive the following system of integral equations for 
the free energy: 

F = -NT J dx [s(x + l/g) + s(x - 1/g)] In \^ 4 - ; ! r 



N T / dx s(x) In 



(8) 



with 



ej (x) = Ts * In 



1 + e *j-i(z)/rl \ 1 + e e j+1 (x)/T 

— Sj^AcOShTTX (9) 



and linin^oo e n (x)/n = H, H being the magnetic field. 
In (|{]), one has A ~ exp(— ir/g), s(x) — [2 cosh ttx] ~ 1 , 
and the convolution product is denoted by *. 

At low temperatures T « A, one has £2(2:) w 
— Acosh7rcc, ei(x) = 0(e _A / T ). In this case the first 
term in the free energy (pf) gives an exponentially small 
contribution corresponding to excitations with a spectral 
gap A and the last term gives F ps —TNq In \/2 showing 
that each central spin contributes ln-\/2 to the ground 
state entropy. A lesson we learn from the exact solution 
is that there are degrees of freedom presumably localized 
on central spins (and their number corresponds exactly to 
a single Majorana mode) which remain decoupled. One 
may conjecture that since the magnetic modes are sepa- 
rated from the ground state by a gap, the soft modes will 
remain soft even when one departs from the integrable 
point. Below we shall give more rigorous arguments to 
support this conjecture. 

The Emery-Kivelson limit. We shall now study 
a U(l) version of the Hamiltonian (^), characterized by 
anisotropic interactions with g — > g« , g± where in the 
limit of strong anisotropy, the low energy degrees of free- 
dom can be identified. This approach has been fruitfully 
applied in quantum impurity problems and gives a sim- 
ple description of the two-channel Kondo model and 
also of the Kondo lattice |ll[]. In particular, for the two- 
channel Kondo problem, Emery-Kivelson identified the 
residual zero point entropy stemming from a unique Ma- 
jorana zero mode. To reformulate the Hamiltonian (||) as 
a fermionic theory similar to Kondo problems, we first in- 
troduce three right- and left-moving Majorana fermions 



2 



Tb, h (b = 1, 2, 3) to use the fermionic representation of 
the SU(2) 2 spin current: J a = — ie abc r;,r c /2 with a simi- 
lar relation for the left-moving current . The next step 
in our solution is to combine the two Majorana fermions 
(ri,/i) and (rj,^) to form a single Dirac spinor (R,L) 
which in turn can be bosonized. Introducing a massless 
bosonic field $ and its dual field Q, one can write: 

fi ^ !1 7r = 7S5 exp[ '^ ( *- e)1 

L = ^ = ^ CXI>H ^ ( * +e)1 ,10 » 

where ao is a short distance cut-off. The anticommuta- 
tion relation between R and L is taken into account by 
the commutator: [G(y), $(2:)] = — i6(y — x). The real 
fermionic zero mode (e) is necessary to establish the cor- 
rect anticommutation relations with the third Majorana 
fermion (r3, £3). The interacting part of the Hamiltonian 
([|) is then given by: 

(Vv^e""/ 4 ^ + e ^ c'" '/•, ,) + fc.c) . (11) 

Following Emery and Kivelson, we absorb the phase fac- 
tor q^V^&j into the spin operators by a unitary trans- 
formation U : 

C^e-^Ei 6 '"?. (12) 



As a result the interaction becomes 




where we have replaced the combinations of spin-1/2 op- 
erators and fermionic zero mode by local fermions (oj) 
using the Jordan- Wigner transformation. We have also 
absorbed a phase factor (— i) J in the definition of Oj , 

At a special point = -k (called the Toulouse point), 
part of the interaction vanishes and the low energy 
physics can be studied by a simple mean-field theory. 
Since the scaling dimension d of the bosonic exponents 
in ( |l3| ) is pretty small (d = 1/4), we shall replace them 
by their averages and try to solve the problem self- 
consistently. Introducing two Majorana fermions (£1,2) 
associated with the complex fermion a: a — (£i+i£ 2 )/ v2, 
we find that the effective action decouples into bosonic 
and fermionic parts in the mean field limit: Smf = 
Sb + Sf with 

S B = J dxdr [i (d^) 2 + Acos (Vtt$)] 
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+ 2*j d rXi - %'+iXj ~ iA 6jXj] (14) 

where \ is a Majorana fermion express in terms of 
the right- and left-moving Majorana fields (^,^3): 
Xj = (h,j + (— l) J ?*3,i)/v2- The mean-field parame- 
ters are defined by: A = — ig± (£2%) / V 7ra o an d A = 
gj_(cos(y / 7r<I > ))/y / 7rao. Solving the self-consistency equa- 
tions, we find that the bosonic field <& becomes massive 
with a gap of the order A ~ g^ \v} ^ (\ / g\aa) . In 
the fermionic sector, the Majorana fermion £ 2 hybridizes 
with the Majorana field \ which in terms of the original 
spins stems from the rails spins. The resulting excitation 
spectrum is reminiscent of the one found in Kondo lat- 
tices Q with a small gap A g ~ A 2 /J|| < A. Finally, 
there is still a singlet localized Majorana fermion £1 which 
decouples from the conduction sea at the Toulouse point. 
This mode gives a zero point entropy of magnitude ^ In 2 
per central spin as found above in the integrable model. 
Away from the Toulouse point (8g\\ = g\\ — n =/= 0), this 
mode will acquire a small dispersion and will contribute 
to the specific heat. Coming back to the original model, 
adding the contribution of the singlet Majorana fermion 
(tq, Iq) which has decoupled from the interaction, the to- 
tal central charge in the long-distance limit is C=l. Let 
us stress that the two Majorana modes, in the singlet 
sector, contributing to the central charge are of different 
nature. The Majorana field (ro,Zo) describes a critical 
Ising model whereas £1 is a strongly localized Majorana 
fermion. In contrast, the triplet sector has a small spec- 
tral gap. As a result of this spectrum, the middle spins 
are disordered and have short-ranged spin correlations. 

The experience gained from the study of Kondo models 
leads us to expect that the preceeding results obtained 
in the Emery-Kivelson limit of the model will extend to 
the isotropic point. To show this, we shall now develop a 
mean-field theory directly at the isotropic limit keeping 
track of the lattice structure more accurately than in the 
previous approach. 

Direct mean-field approach. The Hamiltonian (^) 
is reformulated in terms of a fermionic model on the lat- 
tice: 

n t = iJ\\ ]T x a 0+ iX a 3 " f E ■ (iX2j x X2j] 

j j 

+ [X2j+1 X X2j+l\) ■ (15) 

In the continuum limit, the rail of x a fermions will 
contain left- and right- moving Majorana fermions 
(r ai l a ;a = 1,2,3) of the SU(2) 2 spin current: Xj = 

—l a (x) — (— l) J r a (x) ; with this identification, the Hamil- 
tonian (Q) is then reproduced. The model ( [l5|) corre- 
sponds to local moments {$23+1/2) interacting with a sea 
of three Majorana fermions (Xj)- To describe the mid- 
dle spins we use the Majorana representation for spins 
S=l/2 (see Ref. Q and references therein): 
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s n = -^[%x%], n = 2k+-,k£Z (16) 

where 7^ are local Majorana fermions satisfying the an- 
ticommutation relations {7^,7^} = S nm 8 ab . This rep- 
resentation ( |l6| ) reproduces the spin commutation rela- 
tions and gives the correct value of the Casimir operator: 

S " = 3/4 r-h rtn 

A model very similar to (n3) was analyzed in Rcf. p3[ . 

Following this analysis we substitute (|l^) into (|l^) and 
decouple the interaction with an auxiliary field V living 
on the links connecting rails with middle spins. In the 
mean-field approximation the variables V are considered 
as static with their values determined self-consistently by 
minimizing the free energy. We find that the minimum 
is achieved when a unit cell contains two middle spins 
(see Fig. g) with either U+ = V- = 0, U- = ±V+ = A or 
U- = V+ = 0, U+ = ±V- = A. 



In conclusion, the version of the kagome lattice stud- 
ied in this paper can be reformulated as a fermionic the- 
ory similar to models of Kondo lattices. Its low-energy 
excitations are two spin-singlet Majorana modes with 
different spectra: a critical Ising mode and a strongly 
localized Majorana fermion, whereas the triplet sector 
has a small spectral gap. The physics of the localized 
low-energy mode is similar to physics of the two-channel 
Kondo model. This picture retains some properties of 
the kagome antiferromagnct with very soft singlet modes 
and a gap for magnetic excitations. It is tempting to con- 
jecture that the singlet degrees of freedom and the gapful 
magnetic excitation identified here might be responsible 
for the additional structure seen in the specific heat of 
the kagome magnet at low temperature. 
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FIG. 2. The unit cell in the mean-field approximation. 

All branches of the spectrum are found to be gapful. 
The low energy band is rather flat towards the edge of 
the Brillouin zone in agreement with the Emery-Kivelson 
limit of the model in the triplet sector. However, as was 
demonstrated in Ref . jl3| , the local Z2-degeneracy of the 
ground state generates a local real fermionic zero mode 
70 which is coupled to the three Majorana band fermions 
with the Lagrangian: 

L = Y^^ + HX]X^X% (17) 
3 

This singlet zero mode acquires a small dispersion with 
the bandwidth A g ~ A 2 /J|| <C A. The total amount 
of entropy accumulated in this band is h In 2 per central 
spin and this degeneracy will be slightly lifted and results 
in the coherent dispersion of this mode. This mode corre- 
sponds to the strongly localized Majorana fermion found 
in the Emery-Kivelson limit. This field will manifest it- 
self in the spin-spin correlation functions of the central 
spins since at low energies, they behave like s ~ 70 \. Due 
to the small dispersion of the fermionic mode 70, corre- 
lation functions of middle spins are strongly localized in 
space, but not in time where the characteristic scale is 

(( S (r,j) S (0,j)))^K 1 (A g \r\) (18) 

which for small times is proportional to 1/|t|, like for the 
two-channel Kondo problem. 
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